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Abstract 

We consider the extension of the statistical parton distributions to include 
their transverse momentum dependence and we take into account the effects 
of the Melosh-Wigner rotation for the polarized distributions. With a suit- 
able choice of the fragmentation function, we make predictions for the recent 
semiinclusive Deep Inelastic Scattering data on the cross section and double 
longitudinal-spin asymmetries from JLab. We also give some predictions for 
future experiments on electron-neutron scattering. 
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1 Introduction 



A new set of parton distribution functions (PDF) was constructed in the 
framework of a statistical approach of the nucleon PQ, and let us first re- 
call very briefly, its main characteristic features. For quarks (antiquarks), 
the building blocks are the helicity dependent distributions q±(x) (q±(x)) so 
this allows to describe simultaneously the unpolarized distributions q(x) = 
q+(x) + q-(x) and the helicity distributions Aq(x) = q+(x) — q~(x) (similarly 
for antiquarks). At the initial energy scale taken at Ql = 4GeV 2 , these dis- 
tributions are given by the sum of two terms, a quasi Fermi-Dirac function 
and a helicity independent diffractive contribution, which leads to a universal 
behavior at very low x for all flavors. The flavor asymmetry for the light sea, 
i.e. d(x) > u(x), observed in the data is built in. This is clearly understood 
in terms of the Pauli exclusion principle, based on the fact that the proton 
contains two u quarks and only one d quark. The chiral properties of QCD 
lead to strong relations between q{x) and q(x). For example, it is found 
that the well estalished result Au(x) > implies Au(x) > and similarly 
Ad(x) < leads to Ad(x) < 0. Concerning the gluon, the unpolarized dis- 
tribution G(x, Qq) is given in terms of a quasi Bose-Einstein function, with 
only one free parameter, and for simplicity, one assumes zero gluon polariza- 
tion, i.e. AG(x,Qq) = 0, at the initial energy scale Qq. All unpolarized and 
polarized light quark distributions depend upon eight free parameters, which 
were determined in 2002 (see Ref. pQ), from a next-to-leading order fit of a 
selected set of accurate DIS data. Concerning the strange quarks and anti- 
quarks distributions, the statistical approach has been applied to calculate 
the strange quark asymmetry and the corresponding helicity distributions, 
which were found both negative at all x values More recently, new tests 
against experimental (unpolarized and polarized) data turned out to be very 
satisfactory, in particular in hadronic reactions, as reported in Refs. |3j H] . 
The paper is organized as follows. In the next section we review the con- 
struction of the statistical distributions and we present an improved version 
of the extension to the transverse momentum dependence (TMD). In section 
3, we will consider charged pion production in semiinclusive deep inelastic 
scattering (SIDIS), £N — > £HX, a suitable reaction for testing our TMD 
distributions, more specifically, for the cross section and the longitudinal-spin 
asymmetry, by taking into account the effects of the Melosh-Wigner rotation. 
The results are given and discussed in section 4 and the last section is devoted 
to our concluding remarks. 
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2 The TMD parton distributions 



2.1 The original longitudinal parton distributions 

We now review some of the basic features of the statistical approach, as op- 
pose to the standard polynomial type parametrizations of the PDF, based 
on Regge theory at low x and counting rules at large x. The fermion distri- 
butions are given by the sum of two terms pLJ , a quasi Fermi-Dirac function 
and a helicity independent diffractive contribution equal for all light quarks: 



AX^x b Ax h 
exp[(a; — X^A/x] + 1 exp(x/x) + 1 



Ql) = u '"Tw-n ■ ; + , (!) 



- h( n 2, AjXjrW Ax b 

xq [x, Li ) + x _ h)/s] + i + exp(x/ _ } + 1 , 



at the input energy scale Ql = 4GeV 2 . 

Notice the change of sign of the potentials and helicity for the antiquarks. 
The parameter x plays the role of a universal temperature and X$ q are the 
two thermodynamical potentials of the quark q, with helicity h — ±. The 
eight free parameters^ in Eqs. f|T|2|) were determined at the input scale from 
the comparison with a selected set of very precise unpolarized and polarized 
DIS data pp. They have the following values 

x = 0.09907, b = 0.40962, b = -0.25347, A = 0.08318, (3) 

X+ = 0.46128, Xq u = 0.29766, X~ d = 0.30174, X+ = 0.22775 . (4) 
For the gluons we consider the black-body inspired expression 

xG(x,Ql) = j^jzs r, (5) 

exp(x/x) — 1 

a quasi Bose-Einstein function, with be = 0.90, the only free parameter d, 
since Aq = 20.53 is determined by the momentum sum rule. We also assume 



X A = 1.74938 and A = 1.90801 are fixed by the following normalization conditions 
u — u = 2, d — d = X. 

2 In Ref. pQ we were assuming that, for very small x, xG(x, Qq) has the same behavior 
as xq(x, Qq), so we took be = 1 + b. However this choice leads to a too much rapid rise 
of the gluon distribution, compared to its recent determination from HERA data, which 
requires bij — 0.90. 
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that, at the input energy scale, the polarized gluon, distribution vanishes, so 

xAG(x,Ql) = . (6) 

For the strange quark distributions, the simple choice made in Ref. pQ was 
greatly improved in Ref. |2), but they will not be considered in this paper. 
In Eqs. ffT|2|) the multiplicative factors Xq and (X^) -1 in the numerators of 
the non-diffractive parts of g's and g's distributions, imply a modification of 
the quantum statistical form, we were led to propose in order to agree with 
experimental data. The presence of these multiplicative factors was justified 
in our earlier attempt to generate the TMD [5] , as we will explain now, with 
a considerable improvement. 

2.2 The TMD statistical distributions revisited 

Let us recall that the TMD of the non-diffractive part of a quark distribu- 
tion q of helicity h (first term in Eq. ([1])) was introduced by the following 
multiplicative term 

! (7) 

where Y 0q is the thermodynamical potential associated to the quark trans- 
verse momentum kx and l//i 2 is a Lagrange multiplier, whose value is deter- 
mined by a transverse energy sum rule. 

We notice that this term induces a non factorizable x and kx dependence as 
it is assumed in some other parametrizations. In order to recover the original 
x distributions, the integration of ([7]) over k\ gives the explicit factor 

f°° dk 2 
J Q exp[(k T /xfi 2 - Y£)/x] + 1 

Here Lii denotes the polylogarithm function of order 1, which is known to 
arise from the integral of Fermi-Dirac distributions and is such that 

- LiA-eV) = f°° —, — — r- = In (1 + e y ) . (9) 

Similarly for an antiquark distribution q of helicity — h, according to the rules 
of the statistical approach, one should use the same potential with opposite 
sign, so one gets instead, x/i 2 xln(l + exp[— Y h /x}). In Ref. [5], we made the 
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arbitrary simple choice Y^ q = kX$ q , with k = 1.42, which allows to recover 
the factor X^ q for quarks (see Eq. flTJ), since for large values B of Y^ q /x, 
one has Lii(— exp[Y" ^/a;]) ~ Y^/x, which is proportional to X$ q . However 
this is not suitable to get the factor [AqJ -1 for antiquarks (see Eq. fl2])), 
because Lii(— exp[— Y^/x]) ~ exp[— YqJx] for large values of Y^Jx. In other 
words, the product ln(l + exp[Y" ^/x]) ■ ln(l + exp[— Yq q /x]) does not remain 
independent of Y^ q , as it should. 

Actually, the division by x of the argument of the exponential in the Fermi- 
Dirac expression was not necessary because for the transverse degrees of 
freedom, /i 2 plays the role of the temperature. This feature reflects the fact 
that one should not treat on equal footing longitudinal and transverse degrees 
of freedom. Therefore, for the sake of simplicity, we propose to replace Eq. (J7J) 
by ^ 

exp(fc 2 fxfj.2 _ yfc) + i ' ( 10 ) 

with the corresponding integral over k^, x/i 2 ln(l + exp[F g]). Clearly this 
implies a different normalization for fi 2 and Y^ q . At high kr, Eq. (fTOl) has 
a Gaussian behavior, with a width proportional to fiy/x, at variance with 
the usual factorization assumption of the dependences in x and kx {Sj. The 
product ln(l+exp[Y^]) -ln(l + exp[— Y^ q ]) has its maximum (ln2) 2 for Y$ q = 
and therefore it is stationary around this value. 

So now in order to try to recover the factors X^ q and (A^) -1 in Eqs. f !T|2|) . 
we simply have to choose Y$ q such that ln(l + exp[F g]) is proportional to 
Aq ? and more precisely such that 

\n(l + exp[Y h q }) = kX£ q . (11) 

This way we recover exactly the factors X^ q introduced in Eq. ([1]) for the 
quarks. We take the proportionality factor k = In2/A l d , X^ d being the 
lowest longitudinal potential, so with the value given in (jlj) we get k = 3.05. 
In order to get almost exactly (Aq ? ) _1 for the antiquarks in Eq. (|2]), we also 
assume that the corresponding transverse potential Y^ d is small and fixed to 
the value 0.01. So from Eq. ffTTj) . the values of the other three transverse 
potentials can be obtained and we finally have 

Y+ = 1.122, Y - u = 0.388, Y~ = 0.409, Y+ = 0.010 . (12) 
3 x has a small value according to Eq. (j3|) above. 
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These are different from the values obtained in Ref. [5] and will lead to differ- 
ent predictions for the TMD of the PDF. The non-diffractive contributions 
read now 

m h (r k 2 ) = - f1"Vl 

q { ' T) exp(x - X^/x + 1 exp(k 2 T /x^ - Y h q ) + 1 ' { ' 

^ 1 



where 



xqTix, ki) = ^4 - h , (14) 

V TJ exp(x + X - /l )/x + lexp(A;2/V + F /) + l V ; 

F ( x \ = ±°2 = a*) 

{ } ln(l + ex P y h g )/i 2 kfi 2 ' 1 ' 

Similarly for q we have F(x) = Ax 2b ~ 1 /kfi 2 . After this new determination of 
the transverse potentials, we will see later how we can determine /i 2 , using 
the transverse energy sum rule. 

As noted in Ref. [5], if p z denotes the proton momentum, its energy can be 
approximated by p z + M 2 /2p z , where M is the proton mass. Similarly the 
energy of a massless parton, with transverse momentum kj- is, in the same 
approximation, xp z + k\j1xp z . Therefore all involved parton distributions 
denoted Pi(x, k^) must satisfy the momentum sum rule 

J dx J xpi(x, kfydk^ = 1 , (16) 

and also the transverse energy sum rule 

V / dx I Pi {x, k 2 T ) k idk 2 T = M 2 . (17) 
■ Jo J x 

The contribution of Eq. (113p . for the quarks, to the sum rule Eq. ( TT71) is given 
by: 



^q h (x 1 k 2 T )dxdk 2 T = ^4 dx f 

J x ' J * x 2 (exp^ + l) exp(^-^) + l 

(18) 

and after the change of variable £ = k^/xfi 2 , we get 
* l /j, 2 F(x)dx f°° £d£ 
exp 



x - x °<* + i Jo exp(£ - Yfi) + 1 
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where 

f 1 fi 2 F(x)dx = P qh 
1 Jo exp ^ + i ln(l + expy o y ' 1 ] 

P q h is the number of parton of type q h , and 

Therefore in the limit = the contribution of a parton of type q h is just 
/i 2 P g / l 7T 2 /(121n2), since Li 2 (-1) = -vr 2 /12. 

In a similar way for the contribution to the sum rule Eq. ( TTTj) . from the 
non-diffractive part of the light antiquarks Eq. ( I14p . we get 

? r r r. . . = . & , (22) 



exp 



(k, 



1 J ex P (e+y /) + i 



where 



fl 2 F(x)dx Pgh 



eX pf±JoT + 1 ln(l + exp(-y - ft )) 
is the number of parton of type q h , and 



(23) 



ex P (e+y ^) + i 6 



°9 ' , T it ^^txr-h 



Li 2 (-exp(F - ft )) . (24) 



Finally we turn to the universal diffractive contribution to quarks and 
antiquarks in Eqs. ffTPj) . namely xg D (:r, Qq) = Ax b /[exp(x/x) + 1]. Since 
b < (see Eq. ([3])), the introduction of the k T dependence cannot be done 
similarly to the non diffractive contributions, because in the energy sum rule 
Eq. (jTTl) . it generates a singular behavior when x — > 0. Therefore in order to 
avoid this difficulty, as in Ref. [5], we modify our prescription by taking at 
the input energy scale 

Xq (X ' t) = ln2/i 2 [exp(x/z) + 1] [exp(£; 2 /x 2 fi 2 ) + 1] ' ( } 

whose /ct fall off is stronger, because xfi 2 is now replaced by x 2 /i 2 . Note that 
this is properly normalized to recover xq D (x,Q'fy after integration over k\. 
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We have checked that xq D (x, k?) gives a negligible contribution to Eq. (TT7|) . 
as expected (See Appendix). Concerning the gluon, since it is parametrized 
by a quasi Bose-Einstein function, one has to introduce a non-zero potential 
Yq, in contrast with the QCD equilibrium conditions, to avoid the singular 
behavior of Lii(exp[— Yq/x]), when Yq = 0. The value of Yq is not con- 
strained, but by taking a very small Yq, it does not affect the energy sum 
rule (See Appendix). 

Clearly these regularisation procedures to the diffractive contribution and to 
the gluon are not fully satisfactory, but we will discuss them again in our 
concluding remarks (See Section 5). By summing up all contributions to the 
energy sum rule, one finally gets the value of /z 2 , namely /i 2 = 0.198GeV 2 . 

2.3 The Melosh-Wigner transformation 

So far in all our quark or antiquark TMD distributions (see Eqs. (jl3)14p ). 
the label '"/t"' stands for the helicity along the longitudinal momentum and 
not along the direction of the momentum, as normally defined for a genuine 
helicity. The basic effect of a transverse momentum hr ^ is the Melosh- 
Wigner rotation [7J [5], which mixes the components q in the following way 

q + = cos 2 8 q + + sin 2 8 q~ and q~ = cos 2 8 q~ + sin 2 8 q + , (26) 

where 28 = ArctglKkx/xM), M is the proton mass and k is a dimensionless 
parameter. 

Consequently q = q + + q~ remains unchanged q' = q, whereas we have 

Ag' = (cos 2 6» - sin 2 6») Aq = cos26» Aq = cos Arctg(/tfc T /xM) Aq . (27) 
So we finally get 

Ag' = 1 Ag . (28) 

y/1 + {^k T /xMf 

The effect of the Melosh-Wigner transformation on the double longitudinal- 
spin asymmetry will be discussed in Section 4. 

2.4 The TMD distributions in the relativistic covariant 
approach 

Covariant parton models have been widely discussed in the literature, but 
in some recent papers [9j [10], an analysis based on the requirements of sym- 
metry, for the parton motion in the nucleon rest frame, leads to a different 
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method to generate the TMD of a given x-distribution. By using some in- 
put unpolarized distribution f(x), one can calculate the corresponding TMD 
distribution f(x,k^), by means of its derivative, according to the following 
rule 

= -^</«)/f>. < 29 > 

where the variable £ is defined as £ = x(l + k 2 -.jx 2 M 2 \ M being the proton 
mass. 

This method has been generalized for helicity distributions Af(x) and in this 
case we have for the corresponding TMD distribution Af(x, k\) 

A/( "'^ ) = ^ [3A/(0 + 2 / 1 ^^-4 A/(0] • (30) 

It is interesting to recall that in Ref. [llj, it was demonstrated that for the 
TMD PDF a factorized form fi(x)fz(kj<) is in contradiction with Lorentz 
structure, at least for zero strong interaction coupling g — 0. Using a rather 
different approach, they obtain results identical to the above ones. 
We will show and discuss later the results one obtains from these formulas, 
using as input the x-dependent statistical PDF in Eqs. (1,2). In this approach 
as well as in ours, one gets distributions function of x and k\jx. 



3 Cross section and spin asymmetry of pion 
production in polarized SIDIS 

Following Ref. [12] . we consider the polarized SIDIS, £N — > £H X in the 
simple quark-parton model, with unintegrated parton distributions. Accord- 
ing to the standard notations for DIS variables, £ and £' are, respectively, 
the four-momenta of the initial and the final state leptons, q = £ — £' is the 
exchanged virtual photon momentum, P is the target nucleon momentum, 
Ph is the final hadron momentum, Q 2 = —q 2 , x = Q 2 /2P ■ q, y = P ■ q/P ■ £, 
z = P-P H /P-q, Q 2 = xy(s — M 2 ) and s = (£ + P) 2 . We work in a frame with 
the z-axis along the virtual photon momentum direction and the x-axis in 
the lepton scattering plane, with positive direction chosen along the lepton 
transverse momentum. The produced hadron has transverse momentum 
(For further details see Ref. [T2]). 

Keeping only twist-two contributions and terms up to 0(M/Q), the cross 
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section for SIDIS of longitudinally polarized leptons off a longitudinally po- 
larized target can be written as: 



d 5 a ^ 2a 2 
, , , , 2 = — 2~ {^i + ^ ^lH 2 } , (31) 
ax ay az cl z pt xy z s 

where the arrows indicate the direction of the lepton (— y) and target nucleon 
(<=) polarizations, with respect to the lepton momentum; A, and Sl are the 
magnitudes of the longitudinal beam polarization and the longitudinal target 
polarization, respectively. 

The two terms have the following simple partonic expressions 

Ui(p T ) = X>?/ d 2 k T q(x,k T )^^-D^z,q T ), (32) 

H 2 (p T ) = E e ?/ d 2 k T Aq'(x,k T )7ry 2 ^^D h q (z,q T ), (33) 

where Pt = Qt + zkj- and qx is the intrinsic transverse momentum of the 
hadron H with respect to the fragmenting quark direction. Here s, t and 
u are the Mandelstam variables for the subprocess Iq — > £q. Note that in 
Eq. ( |33i) above, we have used Eq. ( |28i) . which takes into account the effect of 
the Melosh-Wigner rotation. 

The first two contributions, Eqs. (|32|) and (|33|) . give, respectively, the un- 
polarized cross section and the numerator of the double longitudinal-spin 
asymmetry A\ 

d 5 a 2a 2 d 5 a ++ d 5 a + - Aa 2 

rii ~r — ; — ; — « ; — ; — ; — « — — — t~ ^2 ? 



dxdydzd?pT xy 2 s dxdydzd 2 pr dxdydzd?px xy 2 s 

(34) 

where +, — stand for helicity states. So we simply have A\ = 27^2/^1- 

The integrals in Eqs. f!32|l33j) involve the following TMD fragmentation 
function [13] 

D h q (z, q T ) = D h q {z) * exp , (35) 

which is the standard factorized Gaussian model, since we have not yet gen- 
eralized our statistical approach to the TMD fragmentation functions. 
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4 Results and discussion 



In this section we will present all our results on the TMD unpolarized and 
polarized PDF, for light quarks and antiquarks, resulting from the two ap- 
proaches considered above. We will discuss their specific features and the 
difference they lead to, in the calculation of the cross sections and the 
spin asymmetries for SIDIS pion production. All our results are given at 
Q 2 = 2.37GeV 2 , the value corresponding to the CLAS data [HI US], so 
we have performed a backward QCD evolution from our input energy scale 
Q 2 = 4GeV 2 . 

In Fig. [T] we show xu(x, k T , Q 2 ) and xd(x, kx, Q 2 ) as a function of kr for 
different x values, using the TMD statistical PDF constructed in section I2T21 
We have checked that in this x-region, the non-diffractive part of the quark 
distributions largely dominate. Similarly xAu(x, kr, Q 2 ) and xAd(x, kr, Q 2 ) 
are shown in Fig. [2] and we recall that they are independent of the diffractive 
contribution. It is clear that all these kx distributions are close to a Gaussian 
behavior, but with a x-dependent width. The corresponding antiquark PDF 
are shown in Figs. [3] and H] and in this case, we notice a much rapid fall off 
in kx compare to the unpolarized PDF. 

Now if one uses the procedure resulting from the relativistic covariant ap- 
proach described in section 12.4} one obtains different TMD unpolarized and 
polarized PDF as shown in Figs. and [H] for u and d quarks. By comparing 
them with Figs. [1] and [2j we see that their kr fall off is much faster than in 
the previous case. The corresponding antiquark distributions are shown in 
Figs. and |HJ 

Next, we turn to the calculation of the unpolarized cross section and the 
double longitudinal-spin asymmetry for pion production in polarized SIDIS. 
The cross section is directly related to T-L% (see Eq. f )34|) ) and we show in 
Fig. |9j the results, on a proton target, from the two approaches, versus p?, 
for different x- values. The parameter fi 2 D , which enters in the fragmentation 
function, Eq. fl35|) . is a free parameter. In order to get the best description 
of the data of Ref . [T4] . it has been adjusted to the value jJL 2 D = 0.155GeV 2 , 
a value slightly different from the one used in Ref. [5J. The agreement is 
better in the case of the relativistic covariant approach, which has a faster 
•p\ fall off. However it has a very little x-dependence, which is not known ex- 
perimentally at the moment, unfortunately. We predict essentially the same 
result for 7r~ production and also for 7r ± production on a neutron target. 
However if we consider the ratio of the cross sections for the production of a 



11 



7r + from a neutron target over its production from a proton target, our pre- 
diction is shown in Fig. [TUJ The calculation was done in the two approaches 
considered and the results are almost identical. At fixed px, this ratio de- 
creases with x, following almost the trend of the ratio d(x)/u(x) (see Fig. 4 
of Ref. and at fixed x, it is essentially flat over px, because the TMD 
of the pion fragmentation function is flavor independent. This prediction is 
worthwhile to check with future experiments. 

Finally, let us consider the double longitudinal-spin asymmetry A%, defined 
above. The results of the calculation from both approaches, for tt ±, ° on a 
proton target are shown in Fig. [Til with the kinematic cuts corresponding to 
the JLab recent data [T5]. Whereas the relativistic covariant approach leads 
to an asymmetry decreasing with p T , the statistical approach leads to a flat 
dependence in px, in fairly good agreement with the data, and it gives the 
correct normalization. We note that this behavior, which was obtained with 
k = 1.35 in ([28]) . is partly due to the effect of the Melosh-Wigner rotation. 
In both approaches this effect reduces Ai, but it plays an essential role in 
the statistical approach because it compensates the effect of the kx rising 
behavior of Aq. For completeness, in view of future experiments, we have 
also calculated the asymmetry on a neutron target, which are displayed in 
Fig. [121 For the production of tt + , A\ is sensitive to Ad and this is the reason 
for a negative result. In this case, the predictions from the two approaches 
lead again to rather different results, which should be also compared to the 
predictions from Ref. [B|. 

Before closing this discussion we must come back to the effect of the Melosh- 
Wigner rotation. It is clear that the integral over kx of Aq' (see Eq. ([28]) 
is smaller than Aq(x). Therefore to solve this small mismatch, one should 
adjust the potentials, in such a way that Xq — Xq increases slightly, whereas 
Xq q + X^ q remains unchanged, since the Melosh-Wigner rotation does not af- 
fect the unpolarized distribution q. This new improvement will be considered 
more seriously in future work, when we will have access to more precise data 
on the kx dependence of the quark distributions, allowing also a good flavor 
separation between u and d quarks. 

5 Concluding remarks 

An important result of this work is the construction of a new set of TMD 
statistical distributions. This allows us to take into account, in a satisfactory 
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way, the multiplicative factors Xq and (X^) -1 in the numerators of the 
non-diffractive parts of g's and g's distributions. We have introduced some 
thermodynamical potentials Yq , associated to the quark transverse momen- 
tum kxi and related to Xq by the simple relation ln(l + exp[Y" g]) = kX^ q . 
This approach involves a parameter fi 2 , which plays the role of the temper- 
ature for the transverse degrees of freedom and whose value was determined 
by the transverse energy sum rule. The substitution x — > x 2 , we had to 
make in the diffr active part of the quark (antiquark) distributions and in the 
gluon distribution @, to avoid a singularity in the energy sum rule, remains 
an open problem at the moment. We can give the intuitive argument that 
gluons, as photons in a laser, are created mainly in the forward direction. 
The diffractive part, which comes from their conversion into qq pairs, as well 
as the gluons themselves, may thermalize only for the x degree of freedom. 
We have calculated the pr dependence of SIDIS cross sections and double 
longitudinal-spin asymmetries, taking into account the effects of the Melosh- 
Wigner rotation, for 7r ± '° production by using this set of TMD parton dis- 
tributions and another set coming from the relativistic covariant approach. 
These sets lead to different results, which were compared to recent experi- 
mental data. Both sets do not satisfy the usual factorization assumption of 
the dependence in x and kr- We have made some predictions for future ex- 
periments with neutron targets, which will allow futher tests of our results. 
Major progress in our understanding of the TMD PDF will be certainly 
achieved also by measurements at a future electron ion collider [16J. 



A Appendix 

Let us consider the contribution of the diffractive term to the energy sum 
rule. So by using Eq. ([25]) . this contribution reads 

/ %°(x, k 2 )dxdk 2 T = f a ff^' ds f ff- , 
J x J /i 2 x 2 ln2(exp § + 1) J exp (Jk_) + 1 

(A.l) 

4 Let us recall that, unlike the non-diffractive part, they didn't require any multiplica- 
tive factor. 
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and after the change of variable £ = k^/x 2 ^ 2 , we get 

2 f 1 Ax l dx [°° grig fi 2 n 2 P D 
V J expf + lj expg + 1 12 In 2 1 ' ' 

One finds that Pq = 0.0115, so the contribution of the diffractive part to the 
sum rule is negligible. 

Finally let us consider the case of the gluon. In order to recover its in- 
dependence (see Eq. ©), we can write 

XG{X ' ^ = ~ (exp | - l)^ln(l - exp Y G ) exp + Yq) _ \ ' ^ 

because we had to introduce a small Y G , otherwise xG(x, kj) = 0. The 
contribution to the energy sum rule Eq. (fT7|) is 



— G(x, k T )dxdk T — — I — — —dx 



x ' ' Jo V 2 - expF G )(exp | - 1) J exp + y G ) - \ 

(A.4) 

and after the change of variable £ = k^/x 2 ^ 2 , we get 
- ^ / 7 5 TTTn rrr / 77 , „ ; 7 = ^ Pg'^g , (A.5) 



o (expf -l)ln(l- exp Y G ) J exp(£ + Yq) — 1 

0.421 (A.6) 



with 

^ A G x b °dx 

"G 



exp f - 1 



and 



j _ 1 r Li 2 (expF G ) 

G ln(l - exp Yq) J exp(g + K G ) - 1 ln(l - exp Yq) 1 ' ' 

With Y G = 10~ 6 , we find I G = 0.119, so the contribution of the gluon to the 
energy sum rule is negligible. 
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Figure 1: The statistical distributions xu(x, kx, Q 2 ) {left) and xd(x, kx, Q 2 ) 
(right), calculated at Q 2 = 2.37 GeV 2 , versus kx, for different x values: solid 
line x = 0.2, dashed line x = 0.4, dotted line x = 0.6 
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Figure 2: The statistical distributions xAu(x,kT,Q 2 ) (left) and 
xAd(x,k T ,Q 2 ) (right), calculated at Q 2 = 2.37 GeV 2 , versus k T , for dif- 
ferent x values: solid line x = 0.2, dashed line x = 0.4, dotted line x = 0.6 
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Figure 3: The statistical distributions xu(x, kx, Q 2 ) {left) and xd(x, kr, Q 2 ) 
(right), calculated at Q 2 = 2.37 GeV 2 , versus kr, for different x values: solid 
line x = 0.2, dashed line x = 0.4, dotted line x = 0.6 
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Figure 4: The statistical distributions xAu(x,kT,Q 2 ) (left) and 
xAd(x,k T ,Q 2 ) (right), calculated at Q 2 = 2.37 GeV 2 , versus k T , for dif- 
ferent x values: solid line x = 0.2, dashed line x = 0.4, dotted line x = 0.6 
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Figure 5: The relativistic covariant distributions xu(x,kT,Q 2 ) {left) and 
xd(x,kT,Q 2 ) (right), calculated at Q 2 = 2.37 GeV 2 , versus kr, for different 
x values: solid line x = 0.2, dashed line x = 0.4, dotted line x = 0.6 
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Figure 6: The relativistic covariant distributions xAu(x, kT, Q 2 ) {left) and 
xAd(x, hp, Q 2 ) {right), calculated at Q 2 = 2.37 GeV 2 , versus k T , for different 
x values: solid line x = 0.2, dashed line x = 0.4, dotted line x = 0.6 
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Figure 7: The relativistic covariant distributions xu(x,kT,Q 2 ) (left) and 
xd(x,kT,Q 2 ) (right), calculated at Q 2 = 2.37 GeV 2 , versus kr, for different 
x values: solid line x = 0.2, dashed line x = 0.4, dotted line x = 0.6. 
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Figure 8: The relativistic covariant distributions xAu(x, kT, Q 2 ) {left) and 
xAd(x, hp, Q 2 ) (right), calculated at Q 2 = 2.37 GeV 2 , versus k T , for different 
x values: solid line x = 0.2, dashed line x = 0.4, dotted line x = 0.6. 
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Figure 9: The p\ dependence of the term %i at Q 2 = 2.37GeV 2 and z = 0.30 
for 7r + production on a proton target. Comparison of the results of the 
statistical approach (left) and the relativistic covariant distributions (right). 
In both cases the solid lines are for x = 0.20, the dashed lines for x = 0.40 
and the dotted lines for x = 0.60. The data are from Ref. [H] for x = 0.24 
and the error bars are statistical only. 
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Figure 10: The p% dependence of the ratio U^/H^ at Q 2 = 2.37GeV 2 and 
z = 0.30 for 7r + production on a neutron and proton target. Solid line is 
for x = 0.20, dashed line for x = 0.40, dashed dotted line for x = 0.60 and 
dotted line for x = 0.80 . 
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Figure 11: The double longitudinal-spin asymmetry A\ for 7r + (top — left), 
7r~ (top — right) and 7r° (bottom) production on a proton target, versus Pt, 
with the following kinematic cuts corresponding to the JLab data Ref. |15j : 
0.12 < x < 0.48, 0.4 < y < 0.85 and 0.4 < z < 0.7. The data displayed are 
those of Ref. [T5] and the error bars are statistical only. The solid lines are the 
results from the statistical distributions and the dashed dotted lines without 
the inclusion of the Melosh-Wigner rotation. The dashed lines correspond 
to the relativistic covariant distributions and the dotted lines without the 
inclusion of the Melosh-Wigner rotation. 
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Figure 12: The double longitudinal-spin asymmetry A\ for 7r + (left) and 
7r~ (right) production on a neutron target, versus Pt, assuming the same 
kinematic cuts as for the proton target. The solid lines are the results from 
the statistical distributions and the dashed lines correspond to the relativistic 
covariant distributions. 
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